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The exact solution to the Einstein equations that represents a static axially symmetric source 
deformed by an internal quadrupole is considered. By using the Poincare section method we nu- 
merically study the geodesic motion of test particles. For the prolate quadrupolar deformations we 
found chaotic motions contrary to the oblate case where only regular motion is found. We also study 
the metric that represents a rotating black hole deformed by a quadrupolar term. This metric is ob- 
tained as a two soliton solution in the context of Belinsky-Zakharov inverse scattering method. The 
stability of geodesies depends strongly on the relative direction of the spin of the center of attrac- 
tion and the test particle angular momentum. The rotation does not alter the regularity of geodesic 
motions in the oblate case, i.e., the orbits in this case remain regular. We also employ the method 
' of Lyapounov characteristic numbers to examine the stability of orbits evolving around deformed 

nonrotating centers of attraction. The typical time to observe instability of orbits is analyzed. 

^ , I. INTRODUCTION 

(N 

After the visionary work of Poincare [Q and the KAM (Kolmogorov, Arnol'd, and Moser ||^) theory it became well 
^ , established that non-integrability and hence chaos is a general rather than exceptional manifestation in the context 
' of dynamical systems, see for instance Given this ubiquitous fact, an important issue in astronomical modeling 
' is the study in which extent in phase space chaoticity rises in models that are relevant to describe real systems and 
what are its consequences. 

OO The adequate description of the gravitation field of an astrophysical object has been an important subject in both 

relativistic and Newtonian gravity since their origin. The particular case of the gravity associated to an axially 
symmetric body has played a central role in this discussion. Recently, Merrit Q found, from detailed modeling of 
triaxial galaxies, that most of the galaxies must be nearly axisymmetric, either prolate or oblate. In Newtonian theory 
the gravitational potential of an axially symmetric body can be always represented by its usual expansion in terms of 
Legendre polynomials (zonal harmonics). The underlying theory in this case is the usual Newtonian Gravitation that 
O \ for large masses and velocities is known to be less appropriate than Einsteinian General Relativity. In the late case 
• the Newtonian potential is replaced by the spacetime metric and Newton motion equations by geodesies. In General 
^ _ Relativity we have that the solution of the vacuum Einstein equations associated to a static axially symmetric body 
r* ■ has a simple form with only two metric functions |^ and one of them admits an expansion in zonal harmonics. For 
. 5^ ' rotating axially symmetric bodies we have a metric with three functions and two of them obey a sigma-model type 
of partial differential equations for which there are known methods of solution ||^. 
H The change of the particle motion equation and gravitational theory can produce dramatic effects, for instance, test 

- - 1 particles moving in the presence of systems of masses that are integrable in Newtonian theory are chaotic in General 
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Relativity, examples are: the fixed two body problem [[7|Jg], and particles moving in a monopolar center of attraction 
surrounded by a dipolar halo Also gravitational waves, a non existing phenomenon in the Newtonian realm, can 
produce irregular motion of test particles orbiting around a static black hole |l^Jll[]. Another distinctive feature of 
general relativity is the dragging of inertial frames due to mass rotation. This fact is observed, for instance, in the 
impressive differences of the geodesic motion in Schwarzschild and Kerr geometries |]l2| . 

Along this article we shall study the geodesic equations for particles evolving in the space time associated to a 
center of attraction with a quadrupolar deformation. The solution of the Einstein equations representing this center 
of attraction - in the static case - can be found in wherein the rather misleading terminology "distorted black 
hole" was used to refer to such an object. Examples of static centers of attractions with multipolar deformations 
are: a) A true static black hole (or a dense object) surrounded by a distribution of matter like a ring or a small disk 
formed by counter-rotating matter, i.e., built by approximately the same number of particles moving clockwise as 
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counterclockwise. Even though, this interpretation can be seen as a device to have a static stable configuration there is 
observational evidence of disks made of streams of rotating and counter-rotating matter [ p^ . b) An axially symmetric 
static dense object with either polar deformations or polar jets. In the case a) we have oblate deformations. Also the 
polar deformations of the Sun and the inner planets in the solar system are oblate. We have prolate deformations in 
stars with jets and in some galaxy clusters . In the precedent cases, by adding rotation to the central black hole and 
removing the counter-rotating hypothesis we can have stationary centers of attraction with multipolar deformations. 

Geodesic motions for axially symmetric spacetimes representing core-halos system were studied in p6[ | for bounded 
motion and in p7| for unbounded motions. The case of a slowly rotating center of attraction with a dipolar halo was 
considered in |18[ . The geodesic chaos for a disk with a central center of attraction was considered in . A core- halo 
system with NUT (Newman- Unti-Tamburino) charge was also considered . Newtonian and pseudo-Newtonian 
p2| counter parts of some of this systems has been also studied. In a recent paper - within the Newtonian realm - we 
studied chaotic motions of test-particles orbiting around a deformed body modeled by a monopolar and an internal 
quadrupolar term |^^ . 

In this article we dwell in study of geodesic chaos, but now related to internal quadrupole deformations of the 
attraction center. Note that halos are external multipolar contributions, their strength increases with the distance 
contrary to the internal ones that decreases with the distance. The quadrupolar contribution usually represents the 
major deviation to the spherical symmetry. Thus, as a good first approximation, it can model most of the deformed 
sources. 

We shall analyze only bounded motions for specific choices of energy and angular momentum and certain values 
of quadrupolar strength that we believe will cover all the different typical situations. Due to the symmetry of the 
problem, one can reduce the geodesic motion to a dynamical system with two degrees of freedom. For such cases, the 
Poincare section method is the most appropriated tool to study the geodesies general behavior. 

The paper consists of two main parts. In the first one. Sec. II, the exact solution to the Einstein equations that 
represents a static axially symmetric source deformed by an internal quadrupole is considered. By using the Poincare 
section method we numerically study the geodesic motion of test particles. For the prolate quadrupolar case we found 
chaotic motions contrary to the oblate case where only regular motion was found. 

In the second part. Sec. Ill, the rotation of the attraction center is considered. We first study the metric that 
represents a rotating black hole deformed by a quadrupolar term. This metric is obtained as a two soliton solution in 
the context of Belinsky-Zakharov inverse scattering method |Q that generates new solutions from a known one (seed 
solution). As in the precedent section, geodesies were numerically studied using surfaces of section. The consideration 
of different cases leads us to conclude that the black hole rotation considerably alters the stability of the system. We 
obtain that the stability depends strongly on the relative direction of the spin of the center of attraction and the test 
particle angular momentum. We also found that the rotation does not alter the regularity of geodesic motions in the 
oblate case, i.e., the orbits in this case remain regular. We conclude, in Sec. IV, with further considerations on the 
stability of orbits. But, now we employ the method of Lyapunov characteristic numbers following Benettin et al. [p5| . 
We also discuss the typical time to observe instability of orbits and make some final remarks. 



II. SCHWARZSCHILD SOLUTION WITH QUADRUPOLE DEFORMATIONS 

The metric of the spacetime related to the gravitational field of a static axially symmetric source is the one associated 
with the Wcyl line element, 

where (j) Etnd 7 are functions of r and z only. The range of the coordinates r^z^Lp are the usual ones for cylindrical 
coordinates. It is more convenient to use prolate spheroidal coordinates, u and v, that are related to Weyl coordinates 
by 

r2=m2(u2_i)(i_^;2)^ (2) 
z = muv^ 

where m is a constant, that will be associated with the mass of the center of attraction. The coordinate v takes values 
in the interval [—1,1] (it is essentially a cosine) and u runs from 1 to infinity (it is essentially a radial coordinate). 
We shall use units such that c = G = 1. With no lose of generality we shall also choose m = 1. In this new system of 
coordinates, the metric (11]) takes the form , 
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For this line clement the vacuum Einstein equations reduce to, 

[iu'~l)ij,,,U-[il-v^)ij,,U=0, (4) 

_ - 2vi,^,)iu^ - 1)(1 - - m(1 ~ 

_ (2mV>,„ - 7;V;,^)(u^ - 1)(1 - v^)^^^ + v{u^ ~ \fi),u 

Equation (Q) is the usual Laplace equation in prolate coordinates for the metric potential and Eqs. (|^) yield the 
function 7 as a quadrature. The ingrability of 7 (7,111; = l,vu) is guaranteed by Eq. (^. The potential '0 for the 
Schwarzschild solution in prolate coordinates is Q , 

1 , 1 — u 

In this article we shall consider the solution, 



1 1-7/ 

V'=-l0g-— +fc2P2MQ2H, (7) 

2 1 + M 

where Pi and Qi are the second Legendre polynomial and function, 

Pi{x) - (3x2 _ ^)/2, Qi(x) = [P2(:z:)log^ - 3x]/2, (8) 

and ki is a constant that is positive (negative) for oblate (prolate) deformations. Note that the Newtonian limit of 
the potential (g) is ^ = -m/ R+ [2m^k2/lb)P2{cos'd)R-^ . 
From (^ we find the other metric function, 

7 = {4[2((7fc2 - 20A;2 + 4) log(M - 1) + (fcz + 2)^ log(u + 1) - 41og(u2 - w2)(fc2 - 1)^) - 3((27uV - 30^^ - 2lv'^ + 
26)A:2 - 8) log((M + l)/{u - l))k2uv'^ + i{{21u^v'^ - iQu^v^ + iu^ - \2v^ + \%v^)k2 - \'6v^)k2\ - 
3[4((3u2 - 3m - 2)fc2 + 8) - 3(9uV - -u^ _ -^2 ^ !)(„ _ i)(^;2 _ 1) 

log((u + \)l(u - \))k2\{u + 1) log((u + l)/(u - l))fc2}/64. (9) 

The exact solution to Einstein equations given by (0)-(||) was first study by Erez and Rosen |2^, we will comeback 
to this point latter. The general case (Schwarzschild with the whole series of multipoles) was considered by Quevedo 
psf and a simple interpretation in terms of bars was presented by Letelier . This solution has been interpreted as a 
"distorted" black hole in . The study of the associated Newtonian multipoles as well as the relativistic multipoles 
for this solution and other multipolar expansion can be found in The geodesic equations for the metric (||) take 
the form, 

"^-1 „2V; , « r^„2 ,.2x„-2^n -2 ,/.m, (w^ - 1)m 



(9„e2'/' + du{[n^ - 1)(1 - «')e-2'A]) _ ij? [{d^(^ _ ^)] + 



dr2 2e2(-r-'/')(M2 - w2) V " > V y ' -^/J ■ (^2 _ ^2)(^2 _ 1 

(^[9.(7 - V.)] - (^^^3^ j - V ^^^-^^ + (^2„,2)(,2,1) ) ' (10) 

'^^^ _ 1-^^ /fl .2^ , a r^.,2 lA/i ..2^.-2^l^ -2 ^„ ,/.M (u^ - l)v 



(9,e2^ + d,{[u^ - 1)(1 - 7;2)e-2'/']) - (^[9,(7 - V')] 



dr2 2e2(7-V')(u2 _„2) V ''LV /V ^ J; \^L^-V/ r/J (^2 _ „2) („2 _ 1) 

-2.. (^[947 - V.)] + yrr^^ ) ~ u (.2_,2)(,2_i) ) ' (11) 

^^^2^(u.v)^^ L = e-2'^("'''»(u2 - l)(l -t,2)^, (12) 
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where t = s/c = s and the overdots indicate derivative with respect t. E and L are constants of integrations related 
to the test particle energy and the angular momentum, respectively. The metric (0) gives a third constant of motion, 



1 = e^^^^'^H^ - e-2V'(«,")(y2 _ 1)(1 _ ^;2)^2 _ g2(7(«,«)-^(«,«))(^2 _ ^2) f ^ (^3) 

\ — 1 \ — J 

The motion of the test particle is completely determined by the solution of the two second order differential equations 
( p^ ) and (pi]). They define a four dimensional phase space, but the motion constants ( p^ ) and (|l^) tell us that the 
motion is effectively realized in a three dimensional surface. Moreover, these constants allow us to define an effective 
potential like function, 

Thus the motion must be restricted to the region defined by the inequality <i>(w, v) < 0. 

Since the geodesic motion of the test particles is performed in a three dimensional effective phase space an adequate 
tool to study this motion is the Poincare section method. As we mentioned before the sign of the quadrupole parameter 
^2 specifies if the source is deformed in a prolate or in an oblate shape. First we shall study the prolate case, k2 < 0. 

From relation (|l^) we can find the appropriated constants to have a confined motion. We find that only three 
combinations of the constants: E (energy), L (angular momentum) and k2 (quadrupole strength), characterize all the 
possibilities of confinement. In Fig.|l| we present the curve <&(w, v) = for L = 3.32, E — 0.937 and k2 — —5.02. We 
have two bounded regions and two unbounded ones. With the same values to L and E and a small change in the 
quadrupolar constant, k2 — —5.0, we obtain the curve plotted in Fig.|^. The two bounded regions merge in a single 
one. The two escape zones remain unbounded. Finally, in figure 3, we present the curve $(m, u) = for L = 3.8, 
E = 0.9731 and k2 = —1. Now the two zones of unbounded motion merge in a single one and the region of bounded 
motions increases. 

We construct Poincare surface sections for the the three sets of constants indicated above. In Fig.^ we present a 
Poincare section for the two bounded regions of Fig. |^. In the middle bounded region we have a typical picture of 
chaotic motion. However the orbits confined in the right hand side bounded region present regular motion. In Fig.^ 
we show the Poincare section obtained for orbits restricted to the closed surface presented in Fig.|[ It is interesting to 
observe that we have a region of irregular motion in the left hand side of the graphic very similar to the one showed 
in the previous figure and in the right hand side a region of regular motion surrounded by a chaotic one. 

In Fig. ^ we show that the motion in the bounded region of Fig. ^ is regular as in the case of a pure Schwarzschild 
black hole ||l^. These results can be understood by studying the effective potential critical points. We recall that 
a pure black hole (fc2 = 0) with adequate values of the constants E and L has an effective potential with a single 
saddle point. When we add the prolate quadrupolar field ^2 < we have a second saddle point for the value of the 
constants of Figs. ^ and ||. In the third case (Fig. ^) the second saddle point disappears and we end up with the 
same dynamical behavior of the test particles as in the pure Schwarzschild black hole case. 

For the case of oblate quadrupole deformation, i.e., ^2 > 0, we found regions of bounded motion very similar to the 
one presented in Fig. |^. But, we did not find a second as in Figs. ^ and |^. This indicates the absence of a second 
saddle point. We studied surface sections for many different values for E, L and fc2 > 0. We always found regular 
motion. 



$(«, v) = 




III. KERR SOLUTION WITH QUADRUPOLE DEFORMATIONS 

Since the Kerr solution represents a rotating black hole, the addition of an internal multipole term can be used 
to model a rotating star or the core of a galaxy. The black hole rotation produces the pure relativistic effect of 
dragging of inertial frames. Then our main goal, in this section, is to study the influence of the black hole rotation 
on the stability of geodesic motions. Letelier and Vieira ||l^] studied the motion of test particles moving around a 
slowly rotating black hole with a dipolar halo. Now we shall study the case of a central body with arbitrary rotation 
deformed by an internal quadrupole term. 

The metric for a stationary axially symmetric spacetime has the general form. 



gttdt^ +2gt4,dtd4) + g<j,<j,d4>^ - {u^ - w^) 



du2 



dv^ 
1 ^ 1 - t;2 



(15) 
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where gtt, gt<p, gcf>cf> and F are function of the coordinates u, v. 

Behnsky and Zakharov presented a solution generating algorithm for metrics with two independent Killing vectors 
Q. They obtained the Kerr solution by applying the method to the Minkowski spacetime (seed solution). The 
application of this solution generating method to more general seeds was studied in . Using the techniques presented 
in this last article we can easily obtain the metric functions gttiu,v), gt,j,{u,v), gcf,cf,(u,v) and f{u,v) that represent a 
Kerr black hole deformed by multipolar terms. We choose as the seed a metric representing a pure quadrupole. Then 
the the two-soliton solutions give us the nonlinear superposition of a Kerr solution with a quadrupolar field. We find, 

g,, ^ (e^(e2«((2e^^^+2^^ {u' ~ v^) ~ e^" {v^ ~ l)){p + if - e^^^ {u + l){u- l)q')q' 
-e^^'{e^"ip + l)^u + l){u - 1) + e4^i(«2 - l)q^){p + if))/ 
(e2^((2e2^i+2^^ {u + v){u - v) - e"" {v - l)^){p + \)\^^'{u - l)\^)q^ - e^^^ 

{e'^ip + + 1)2 + e^^i {v + l)2q2)(p + 1)2)^ (16) 



gt^ = {^2e"{e''"{{2e^^^+^^^{u^-v^) 

-e^"{v^ - l))(p+ 1)2 - e^^^iu + l){u - l)q^)q^ - 
e^^^e^"{p + l)2(u + l)(ti - 1) + e^^i {v^ - l)q^){p + 1)2 + 
(g2Fi (g4f.(p ^ ^)2(,^ + l)(w + 1) + e^"{u - l){v - l)q^){u -v)~ 
e2^^(e2^(p + 1)2(m + l)(t; - 1) + e^^i {u - l){v + l)q^){u + v)){p + l)p)q)/ 

(e2ff ((2e2^i+2J=^^ (^2 _ ^2) _ g2H(^ _ ^)2)(p ^ ^)2 _ ^AF, _ ^^2^2)^2 _ 

e^^'{e^"{p + l)\u + 1)2 + e^^i {v + l)\^){p + 1)2), (17) 



9<l^<l^ = — (18) 

9tt, 

= -(exp[((((4(2 log(-u + 1) + ^lu'^v^ - QQu^v^ + 9^2 - 36i;'^ + 48^2 - 8 log(7i + w) + 

14l0g(w - 1) - 81og(M - V)) + 9(9u2w2 -u^ -v"^ + l)(w2 - l)(w2 - 1) log((u + l)/{u - 1))2 

-12(27uV - 30u^t;2 + - 2 - (21^2 - 5)(„2 _ i)^) \og(^(u + l)/{u - l)))k2 - 
16((3m2 _ 1) \og{{u + l)/(u - 1)) - 6w)(3w2 - I))fc2)/128)](e2^((2e2^i+2^2(u2 - v"^) - 



e 



^"{v - l)2)(p + 1) + e^^i {p - l){u - l)2)(p - 1) + e4^^(e2^(p + l){u + 1)2 - 



where 



{p -1){V+ l)2)(p + l)))/(4e2^^+2i.^+2i^(^2 „ ^2y )^ (19) 

Fi = (-(2(2(log(u + 1) - 3^2 - 31og(u - 1) + 2 log(u - w)) + 3(3w - l)(w + l)u) 

+3{{v + 3 + 2(t; + l)u){v - 1) - (3w - l){v + l)u^) \og{{u + l)/{u - l)))fc2)/16, (20) 

Fa = (-(2(3((3u- l)(i; + l)?i + 2w2) -41og(u + w) +41og(w- 1)) - 

3((3w - l)u -{v + l)){u + l){v + 1) log((ii + l)/{u - l)))fc2)/16 (21) 

H = (((3^2 - 1) log((u + l)/{u - 1)) - 6u){3v^ - l)fc2)/8fc2)/8. (22) 

The quadrupole strength is fca , g is the source angular momentum per the square of the mass and p is defined by the 
relation p^ + q^ = 1. The metric presented above is a particular case of the general solution that represent a Kerr 
metric embedded in a field of multipoles, see for instance |Q and |^^, we will back to this point latter. 

When one performs the limit fca — > in the solution presented above one obtains the Kerr metric in Boyer Lindquist 
coordinates, r and that are related to the prolate spheroidal coordinates, u and v hy u = (r — m)/a and v = cos 9. 
The constant p and q are related to the Boyer Lindquist constants hy p — a/m, q — a/m, and m"^ = a'^ + a} . m 
is the monopole mass, a is only an auxiliary constant and a is interpreted as the black hole angular momentum per 
unit of mass measured by a very distant observer. 

As in the precedent case, the geodesic equations have two constants of motion, L and F, that obey the relation. 

gut + fft<ji(/> = g^^cj) + gt^i = L. (23) 
The evolution equations for u and v are 
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CpU IT — \ 

(Pv 1 — 

" 2r(u2 - i;2) 

2r 



duTiu^ - 1) 



(7i2 — ?j2) 



2r (?i2 -'y2)(u2 „ 1) 
(u2 - l)u 



2r(v2-l) (^^2 _„2)(^2 _ 1) 



-2uw 



a„r(w2 - 1 



(u2 — u2) 



a„r _ (n2 - 1)^; 

^ ~ (w2 -i,2)(^2 _ 1) 
(u2 - 



2r(u2_l) (u2-z;2)(u2_l) 



(24) 



(25) 



As in the precedent case we have two second order evolution equations (24) and (^5|) for the variables u and w, the 
metric and the constants ( p3[ ) give a new constant relating these two variables, 

(26) 



1 = g"£:2 ^2g*'t'EL + g'l"l'E^ - {u^ - v^) 



1*2 — 1 1 — ti2 



In other words, despite algebraic complications, we have exactly the same dynamical situation as before. The particles 
move in an effective three dimensional space. Thus we shall analyze the motion of test particles moving in the 
gravitational field of a rotating prolate deformed body using Poincare sections as in the non-rotating case. 

Since the main new ingredient in the new system is the rotation of the source we shall keep the angular momentum 
L, the energy E and the quadrupole strength ki fixed and we shall consider test particles moving with angular 
momentum parallel to the spin source (co-rotation) and with angular momentum anti-parallel to the spin source 
(counter-rotation) . 

In Fig.0 we present the region of bounded motions for counter-rotating orbits, qL < 0. We take E = 0.93715, \L\ ^ 
3.322, and ^2 — —5.08, and for the rotation parameter, q = 0.002. We notice a situation similar to the one presented 
in Fig.||. The bounded regions for the co-rotation case, > is shown in FigJ|. We see two relatively small and 
distant closed surfaces. The Poincare section for the counter-rotating case is presented in Fig.^. Chaotic motion can 
be perceived in the left hand side of the graphic and in the external part of the right hand side as in Fig.||. In Fig.p^ 
we present the surface section for the same parameters as above, but qL > (co-rotation). We do not find chaotic 
motion in this case. 

We were not able to obtain bounded motion for both large quadrupole strength and large rotation parameter. We 
studied bounded systems with large rotation speed (of the order of 0.1) but with quadrupole strength always less 
than unit. In these cases the study of Poincare sections leads to regular geodesic motion for co-rotation, as well as, 
counter-rotation. We also found that the confinement region may suffer an appreciable change in size and shape. 



IV. DISCUSSION 



As we said before, the exact solution to Einstein equations presented in the two precedent sections are not new and 
different versions of them have already appear in the literature. We have presented them in this work for two reasons: 
a) For easy reference, and b) Mainly, because for numerical analysis we need a faultless solution. The one presented 
here were derived using algebraic computation and checked using the full vacuum Einstein equations in each case. 

Besides the Poincare section, we have another technique to quantify geodesic instability: the Lyapunov characteristic 
number. A/", defined as, 

A/" = hm 

So 
r — > oo 

where a-nd S are the deviation of two nearby orbits at times and t respectively. Using the technique suggested 
by Benettin et al. p^ ] - who studied numerical problems in the computation of Lyapunov exponents and Kolmogorov 
entropy - one can get the largest A/". 

Using the constants that define the bound region presented in Fig. |l| we obtain J\f = (9.0 ± 1.0) x 10^^^, where 
the maximum Af was obtained for u — 2.6, v = and py = 0. In the right region plotted in Fig. ^ we get always 
TV < 5 X 10^^^. This shows, as expected, that the Lyapunov coefficient in the regions with chaotic motion does not 
vanish while in the region with irregular motion it could vanish. Nevertheless, the value of the coefficient for chaotic 
motions is very small, it means that the dispersion of the orbits is slow if compared with local fluctuations in the 
mean potential. 



log(<5/<5o) 



(27) 
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We found chaotic geodesic motions for the system black hole plus internal quadrupole for a very small range of 
parameters, specifically when there is a second bounded region. When one considers a rotating source one just re- 
obtained the same behavior studied in for rotating centers of attraction with halos. The orbits of counter-rotating 
particles are more unstable that the orbits for corotating particles. In reference ||l^ the case of slow rotation was 
considered. 

We also studied some bounded chaotic motion for large rotation speed {q > 0.1) in a Kerr halo system from an 
exact solution that represents an external dipole plus a Kerr black hole. Using the techniques presented in this paper 
we conclude that the irregularity introduced by external multipole terms are much larger than the ones introduced 
by internal multipole terms. We found also that chaotic motions for large rotation speed are more frequently, i.e., we 
have irregular motion for a larger range of the constants E and L. 

It is not easy to predict the role that chaos could have in measurable characteristics of galaxies. Let us choose a 
parameter T as the time to characterize the dynamics of the system e.g., to draw the sections of an invariant torus 
in a regular motion or a chaotic region in the Poincare section for irregular motion. We shall consider that this is the 
minimum time to have observable effects. 

Rearranging the units to observable parameters, we obtain T in years from the expression 

T X Af X 1.6 X lO^^^yr, 

where N is the number of steps in the simulations and M is the mass of the central black hole (in solar masses) . The 
step N varies for different systems. It is about lO'' for a black hole plus halo system and about 10^^ for a black hole 
plus internal quadrupole system. For a typical galactic bulge we have M = lO^^M©. Then for black hole plus halo 
system we have T = lOOOyr which is a very small value when compared with galactic ages and for the black hole 
plus internal quadrupole system, T — lO^^yr which is a very large value compared with the Universe age (lO^^yr). 
Consequently, chaotic relativistic effect may show up in the formation of structures in a black hole plus external 
halo system, work along this line will soon be reported. The relativistic effect due to the rotation of the source 
may be important. The internal deformations do not have a significant contribution in this case. Another possible 
observational manifestation of chaos was studied in [^3| 
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FIG. 1. Boundary contour ($ = 0) for a static black hole 
+ quadrupole system L — 3.32, E — 0.937 and = —5.02. 
There are two escape zones in the left hand side of the picture 
that corresponds to small values of u and a couple of closed 
zones of bounded motion. 



•f f-J. 



FIG. 4. Poincare Section for the values defined in Fig.gj. 
We see chaotic behavior in orbits confined in the first zone 
of bounded motions. But, the motion in the second zone is 
regular 




FIG. 2. Boundary of contour for the same values for L and 
E but k2 = —5.0. The two escape zones remain. However, 
the two closed zones merge in only one. 




3.5 4.5 
U 



FIG. 5. Poincare Section for the values defined in Fig.g. 
We see chaotic motion in the left hand side of the figure and 
in a small external region of the right hand side. 




FIG. 3. Boundary contour for L = 3.8 and E = 0.973 but 
k2 — —1.0. No vestige of the second bounded region is left. 




FIG. 6. Poincare Section for the values defined in Fig 
We have regular motion. 
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FIG. 7. Boundary contour of the Kerr+Quadrupole sys- 
tem for L = -3.322, E = 0.93715, k2 = -5.08 and q = 0.002. 



FIG. 9. Poincare Section for the same value of the param- 
eters of Fig. 1^ We have chaotic motion mainly in the left 
hand side of the picture 




FIG. 8. Boundary contour of the Kerr-|-Quadrupole system 
for L = 3.322, E = 0.93715, = -5.08 and q = 0.002. 
The confinement region is separated in two. They are much 
smaller than in the precedent figure 
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FIG. 10. Poincare Section for the same value of the pa- 
rameters of Fig. ^. We do not see irregular motion in both 
regions 
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